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! Abstract 

We use a new technique to search for faint red stars in the Hubble Deep Field 
(HDF) imaged by the Wide Field Camera (WFC2) on the Hubble Space Telescope. 
We construct a densely sampled stellar light profile from a set of undersampled 
images of bright stars in a low-latitude field. Comparison of this stellar light profile 
to densely sampled profiles of individual objects in the HDF constructed from 
multiple undersampled dithers allows us to distinguish unambiguously between 
stars and galaxies to / = 26.3. We find no stars with V — I > 1.8 in the outer 90% 
of the volume probed. This result places strong and general constraints on the / 
band luminosity of the constituents of the Galactic dark halo: 



X : irn 5, /,O.5M 



Ml > 15.9 + - log^/^j (V-I> 1.8), 

where M is the mass of the objects and / is their density as a fraction of the 
local halo density, taken to be po — 9 x 10 -3 Mq pc~ 3 . If the halo is made of white 
dwarfs, this limit implies that these objects have My 18.4 and V — I <: 2.5. That 
is, they are ^ 2 magnitudes fainter than the end of the disk white dwarf sequence. 
Faint red dwarfs account for < 1% of the Galactic dark halo for Mj < 14, and 
< 6% for Mi < 15 at the 95% confidence level. The density of inter-galactic Local 
Group stars is at least a factor 3000 smaller than the density of local Population 
II stars. 

Subject Headings: dark matter - gravitational lensing - stars: low mass, brown 
dwarfs, luminosity function 
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1. Introduction 



The detection of ~ 7 candidate microlensing events by the MACHO collabora- 
tion (Alcock et al. 1993, 1995, 1996a,b) during two years of observations toward the 
Large Magellanic Cloud (LMC) implies that a significant fraction (0.2 ^ / ^ 1) of 
the dark halo of the Galaxy may be in the form of massive compact halo objects 
(MACHOs). The mean observed time scale of the events is (t e ) ~ 37 days, where 
t e is the Einstein radius crossing time. This time scale is substantially longer 
than would be expected if the MACHOs have sub-stellar masses, M < 0.1 M , 
and if their velocity and spatial distributions were in accord with a standard halo 
model. Within the context of such a model, the best estimate for the mean mass 
is (M) ~ 0.4 Mq, with a range 0.1 & (M) /M© ~ 1 at the 3 a level. The EROS 
collaboration (Aubourg et al. 1993, 1995; Ansari et al. 1996), has detected two 
candidate events from a somewhat smaller data set and with a somewhat shorter 
mean time scale (t e ) = 26 days. The EROS data are consistent with the results 
inferred from the MACHO data. In particular, the EROS experiment was less 
sensitive than MACHO to the relatively long events of the type that MACHO de- 
tected. These long duration events contribute heavily to both the high estimate of 
the total MACHO density and the high inferred mean mass of the detected objects. 

If there is a substantial population of MACHOs, then it is possible that some 
of these objects may be luminous enough and close enough to be directly detected 
in optical light. Red dwarfs were a plausible candidate of this type. Bahcall et al. 
(1994, hereafter Paper I) searched for such objects in deep images taken by the 
Wide Field Camera (WFC2) on the Hubble Space Telescope (HST) to a limiting 
magnitude of / = 25.3. From the null results of that search, they concluded that 
faint red dwarfs (Mj < 14) make up < 6% of the Galactic halo. With additional 
plausible assumptions, Graff & Freese (1996) have derived even stronger constraints 
using the same observational results. 

The new results of MACHO (Alcock et al. 1996b) may point in a different 
direction: white dwarfs (WDs). Although several arguments appear to constrain 
the WD contribution to the halo to be < 10% (Chariot k Silk 1995) or < 25% 
(Adams & Laughlin 1996), the fact that MACHO may be detecting objects in the 
WD mass range means that this candidate must be taken seriously. More generally, 
the nature of the detected objects, whatever they are, can be probed or constrained 
by searching for intrinsically faint stars. 

The Hubble Deep Field (HDF) (Williams et al. 1996) taken with WFC2 on HST 
provides a unique window on the universe (Bahcall, Guhathakurta, & Schneider 
1990; Abraham et al. 1996; Colley & Rhoads 1996). The extreme depth of the HDF, 
which has an equivalent exposure time ~ 10 greater than the field analyzed in Paper 
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I, provides an unprecedented opportunity to find faint stellar objects. The principal 
advantage of going deep is that it allows one to search for faint stars in regions 
of the color-magnitude diagram (CMD) which are virtually devoid of the stars 
that populate the standard Galactic components. The lack of ordinary disk and 
spheroid stars at very faint magnitudes is a result of the finiteness of the Galaxy. 
Thus, by restricting attention to objects within 1.67 mag of the magnitude limit, 
one probes 90% of the available volume for intrinsically faint stars while eliminating 
nearly all of the stars from previously well-studied Galactic populations. Of course, 
by going deep one also increases the total volume probed (for candidate objects of 
fixed luminosity) but this advantage is secondary: a survey of 16 fields each 2 mag 
less deep would have the same volume and would require only 40% of the telescope 
time. However, these fields would contain of order 30 dwarf and sub-dwarf Galactic 
stars within 1.67 mag of the magnitude limit, where we have made the estimate 
based on actual star counts in the HST Large Area Multi-Color Survey ("Groth 
Strip", / = 96°, b = 60°). Hence, the HDF provides a truly unique opportunity to 
search for intrinsically faint stars. 

In § 2, we describe our technique for discriminating stars from extended objects 
to a magnitude limit of / = 26.3. In § 3, we discuss our selection criteria and report 
that no stars are detected. In § 4, we show that the lack of detections means that 
if the Galactic halo is composed of WDs of mass M = 0.5 M , these must have 
Mj > 15.9 (or My > 18.4). For red dwarfs and brown dwarfs of mass 0.08 M , 
the limit is stronger, Mj > 17.2. 

2. Finding Stars in the HDF 

Objects can be detected in HDF down to / ~ 28. The overwhelming majority 
of these very faint objects are galaxies. Hence, for purposes of studying galaxies, 
no star/galaxy discrimination is required. For star counts however, it is essential 
to distinguish unambiguously between the handful of stars being detected and the 
thousands of background galaxies. Because distinguishing stellar from extended 
profiles requires ~ 5 times more photons than just detection, the magnitude limit 
for the star counts will be much brighter than the detection limit, as we quantify 
below. 

In order to classify objects on WFC2 images, we developed and tested in Paper I 
an effective procedure for separating stars and galaxies based on the radial profiles 
of the objects. The empirical stellar radial profile was determined from a large 
number of stars that appear on a WFC2 exposure at low galactic latitude. Many 
stars (falling at random places relative to the pixel grid) were required because 
the WFC2 poorly samples the point spread function (PSF). (The full width half 
maximum of the PSF is approximately 1.2 pixels. See Fig. 1, below.) 
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The exposures of the HDF were taken in groups at several different positions 
on the sky (these groupings are termed "dithers"), each offset slightly relative to 
the HDF field center. This procedure improves the smoothness of the flatfielding 
corrections and helps improve the spatial resolution. The HDF team released 
special "drizzled" images in four filters, F300, F450, F606, and F814. These images 
were produced using a "drizzling" technique (Fruchter & Hook 1996), which takes 
into account the shifts and rotations between the individual dither positions and 
geometric distortions at the image plane, while preserving the flux. 

Since drizzling involves the rebinning of undersampled data, the radial profiles 
of sources in the drizzled image are slightly different from those in the original 
images. Hence, to find the stars on HDF we developed a procedure that combines 
the depth of the drizzled image with the well-understood PSF properties of the 
raw data. 

We obtained the flatfielded and bias-subtracted frames from the Space Tele- 
scope Science Institute. The observing log indicates that some of the exposures 
were affected by significant scattered light from Earth, which appears as a charac- 
teristic X-pattern in the images. The HDF team ignored these frames in construct- 
ing their first-release drizzled images, and we followed their example in this respect. 
(By exposure time, 6% of the F606 and 25% of the F814 frames suffer from signif- 
icant scattered light). Because of the pointing and tracking accuracy of HST, the 
images taken at each dither position are precisely aligned. This alignment makes 
stacking and cleaning (i.e. removing hot pixels and cosmic rays from) the images 
in each dither straightforward, and we used the same procedure as described in 
Paper I to create cleaned images at each dither position. The total exposure time 
in F814 was 92,200 seconds in 39 separate exposures, taken at 8 separate dither 
positions. In F606 the total exposure time was 110,050 seconds in 105 separate 
exposures taken at 11 separate dither positions. 

We searched the drizzled F814 band image for sources down to a magnitude 
limit of / = 26.5, beyond which the images on the dithered images become too 
noisy to classify reliably. We required that the ratio of peak to total flux be at 
least half as large as the value found for known stars. This initial selection yielded 
629 objects in all. Each object was then located on the 19 dither images (8 in F814 
and 11 in F606). The radial profile of the object was determined independently on 
each of the 19 dither images, and these were plotted against the empirical stellar 
profile for the filters F814 and F606 (as determined in Paper I). 

The key to star/galaxy separation is that we can locate the center of a stellar 
image to an accuracy of better than 0.1 pixels, using the symmetry of the stellar 
image about its center. Thus, for any given star and on any given dither, the 
radial profile of the image is sampled in many neighboring pixels; the separations 
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Figure 1. (a) A densely sampled light profile of a test object (circles) is compared to the template profile 

of a star (solid line) that is constructed from a WFC2 image of stars in a low-latitude field (see Paper I). 

Each circle represents one pixel in one of the 8 WFC2 dither images of HDF taken with the F814 filter. 

The distance (in pixels) from the center of the star to the center of the pixel is plotted against the fraction 

of all star light from the dither image falling on that pixel. Note that even though the individual dither 

images are highly undcrsamplcd, the combined profile is well sampled to well within 1 pixel. The close 

agreement with the template indicates that the object is stellar, (b) Same as (a) except that only one 

of the eight dither images is shown, (c) Profile constructed from the single drizzled image. Vertical and 

horizontal axes have been adjusted to account for the difference in pixel size. Note that the profile is nearly 

identical to the template profile except within 0.6 pixels of the center. 
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correspond to different distances between the stellar center and pixel centers. We 
construct diagrams that give the light intensity of images in pixels whose centers 
are located at different distances from the center of the stellar light. We then 
superimpose the diagrams from the separate dithers and thereby create a densely 
sampled profile. 

Figure 1(a) shows an object of moderate signal to noise, which we classify as 
a star (this object appears on Chip 4 of the drizzled images at (x,y) = (350,598) 
and has / = 23.71, V — I — 1.32). The densely sampled radial profile from the 
8 independently centered dither images (shown as circles) matches the empirical 
stellar profile (shown as a line) very well. The radial profile is well sampled because 
the center of the object falls on the various dither images over a good range of 
spatial positions relative to the pixel grid. Figure 1(b) shows the profile of this 
star from a single dither. Note that for only three pixels does the pixel center fall 
within one pixel of the center of the star. However, these three points are sufficient 
to distinguish stars from galaxies provided that the galaxies are extended by ^ 1 
pixels, i.e. > O'.'l. We used plots of this type in Paper I and in Gould, Bahcall, 
& Flynn (1996) to identify stars in 22 WFC2 fields. In Figure 1(c), we show the 
profile of the same object as seen in the drizzled image, where we take into account 
the different pixel size (drizzled pixels are a factor of 0.402 smaller). The drizzled 
profile begins to deviate from the WFC2 profile at about 0.6 pixels. 

We compared the densely sampled profile with the drizzled profile for several 
dozen stellar and nearly stellar objects. We found that for objects / ^ 25.5, 
stars could be easily distinguished from galaxies using either profile. However, for 
/ 25.5 we found several objects whose densely sampled profiles are clearly non- 
stellar, but which could not be distinguished from stars using the drizzled profile. 
Because of the difficulties in classifying objects on the drizzled image, we identify 
candidate objects on the drizzled image but do the classification on the 19 dither 
images. 

Diagrams similar to Figure 1(a) were made for each of the 629 objects found on 
the drizzled F814 frame, and CF and AG independently classified them as either 
"star", "galaxy", "quasi-stellar", or "?". The designation "quasi-stellar" means 
that the object is clearly not a star, but deviates from a stellar profile only within 
1 pixel. 

Occasionally, one or more of the 19 images could not be used for accidental 
reasons (like a bad pixel) and those cases were handled separately. Photometry was 
carried out after the classification, using the same aperture size and calibration from 
F606 and F814 to /, V — I as described in Paper I. Comparison of the classifications 
by CF and AG showed that we could classify stars and galaxies with confidence 
to / = 26.2. We found a total of 17 stellar objects to this limit. In addition, we 
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found 1 blue stellar object beyond the magnitude limit (I = 26.5, V — I = 0). 
All stellar objects satisfying the / band magnitude limit were easily detected and 
photometered in V. 

The photometric errors, as determined from the scatter of the photometric 
measurements made from different dithers, is ~ 0.08 mag at the magnitude limit 
and ~ 0.03 mag at I — 25. The uncertainties induced by photometric errors are 
therefore several orders of magnitude smaller than the Poisson uncertainties. 

3. Selection Criteria 

Figure 2 is a CMD of the stars detected in the HDF together with the adopted 
selection criteria, which are described below. 

First, we restrict attention to red stars, V — I > 1.8, which includes M dwarfs, 
and is expected to include brown dwarfs and old WDs. The local density of old 
disk WDs has been measured by proper-motion studies to My = 17 (corresponding 
to V - I ~ 1.8) (Liebert, Dahn, & Monet 1988). Halo WDs, which are expected 
to be fainter than this limit, would not have shown up in this study because their 
proper motions ji would generally exceed the selection limit of ji < 2. // 5yr~ 1 (or 
perhaps somewhat smaller). 

Second, we set the magnitude limit at J ma x = 26.30. As discussed in the 
previous section, we find that our star/galaxy separation is reliable for I < 26.2. 
The brightest object above this limit that lies in the adopted color range and is not 
obviously a galaxy has / = 23.39. We therefore have detected all stars to / = 26.3. 

Third, we establish a lower magnitude limit J m i n so as to exclude ordinary 
Galactic stars. As discussed in the introduction, the HDF is so deep that one 
expects very few ordinary stars near the magnitude limit. By setting 7 m i n = 
7 max — 1.67 = 24.63, we therefore exclude the regions of the CMD that are heavily 
populated with previously well-studied stars, while preserving 90% of the total vol- 
ume. Figure 2 shows that, with the above-described criteria, the color-magnitude 
diagram is devoid of stars in the region, 24.63 < I < 26.30, V — I > 1.8 

Finally, we note that for V—I ^ 0.7, there are ~ 60 compact non-stellar objects 
that lie near the magnitude limit. Since some of these objects deviate from point 
sources only in the inner fraction of a pixel (< 0."1), one must wonder whether 
there are not other objects of this class which appear perfectly stellar. Indeed, 
we find three such faint stellar objects with A-star colors and / ~ 26 (see Fig. 
2), and are currently investigating the nature of these extremely compact objects 
as a whole. The presence of this class of object complicates the search for faint 
blue stars. However, since they are well blueward of the "halo zone" , they have no 
impact on the principal results of the present paper. 
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Figure 2. Color-magnitude diagram (V — J vs 7) of stellar objects detected in HDF. The errors are not 
shown since they would be of order or smaller than the size of the points. The selection criteria dcliminating 
the "halo zone" are shown as a box at the lower right. These are I < 26.30 (magnitude limit for star/galaxy 
discrimination), I > 24.83 (eliminate ordinary Galactic stars while keeping 90% of the volume probed), 
and V — I > 1.8 (color region of interest for dark halo candidates). 



4. Limits 

Since we detect no objects with V — I > 1.8, we can rule out at the 95% 
confidence level any model that predicts 3.0 or more detections in this regime. In 
particular, since the HDF goes 1 magnitude deeper than the field used in Paper I, 
we can immediately extend the results derived there. For red dwarfs brighter than 
Mj = 15 (the faintest red dwarf ever seen, Monet et al. 1992) the halo fraction 
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must be / < 6%. For Mj < 14, the fraction is ^ 1%. 

To further interpret these results, consider a class of objects each with mass 
M that comprise a fraction / of the dark halo, and so have a local number density 
n = fpo/M where po = 9x 1CT 3 M pc -3 is taken as the local halo density (Bahcall, 
Schmidt, & Soneira 1983). Now suppose that these objects are detectable out to a 
distance d within the HDF with angular area Q = 4.4arcmin 2 . Then the number 
of expected detections iVexp is 

where we have assumed initially that the halo density is uniform over the region 
probed. In fact, equation (4.1) shows that one can begin to place limits on objects 
that can be detected to d ~ 1-2 kpc, where the exact distance limit d depends only 
weakly on the details of the model. For simplicity, we adopt a halo with density 
oc r~ 2 which has a local density p^ = 0.86po & t a distance of 1.5 kpc in the direction 
of the HDF (/ = 126°,6 = 55°). Since the magnitude limit is I < J max = 26.30 and 
since 10% of the volume is excluded by selecting / > I m [ n = 24.63, we find a limit 

Mi > 15.9 + | log (/^T^) (V-I>1.8). (4.2) 

That is, for objects brighter than this limit there are more than 3 expected detec- 
tions, contrary to the observations. 

For WDs of mass M = 0.5 M Q and V - I > 1.8, these limits imply Mj > 
15.9 + 1.67 log(/). To interpret this limit as a limit on the halo fraction, we assume 
that the observed linear color-magnitude relation for the red end of the disk WDs 
(Monet et al. 1992), Mj = 12.9 + 1.2(V - I), can be extended to the fainter halo 
WDs. We then find limits on the halo fraction / of 

/<0.31xlO°- 72 ^- J )- L8 ]. (4.3) 

Hence, under these assumptions, a full WD standard halo is ruled out for V — / < 
2.5 (My < 18.4) and a 33% WD standard halo is ruled out for V - I = 1.8 
(My = 16.9), both at the 95% confidence level. 

We now compare the sensitivity of our results with the measurements of and 
constraints on luminous halo objects obtained with a variety of techniques as re- 
viewed by Mould (1996). Generally, if one is interested in objects that are at 
least as bright as the end of the locally observed sub-dwarf sequence (My ~ 15, 
V — I ~ 3), then there are several probes that are at least as sensitive than the one 
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presented here. For example, Dahn et al. (1995) use proper-motion selected stars 
to construct a luminosity function. Based on this model, we expect to find ^ 1 
star with V — I > 1.8. Similarly, the deep color-selected ground-based survey of 
Boeshaar et al. (1994) can detect such objects over ~ 3 times the volume probed 
by HDF (although such surveys are subject to significant contamination by faint 
galaxies). The real strength of the HDF is its sensitivity to intrinsically faint ob- 
jects, near the limit (4.2). These would have avoided detection in proper-motion 
surveys. For red objects (V — I k, 3.5), the volume probed by HDF is larger than 
for any ground-based photometric surveys and, of course, HDF is free of galaxy 
contamination. 

Finally, we note that the HDF star counts can be used to constrain the density 
of Local Group stars. Local Group giants and sub-giants Mj ^ 1.5, 0.6 ^ V — I ^ 
1.5 could be seen to a distance d ~ 0.9 Mpc. In the outer 90% of this volume, 
there are no more than 2 such stars observed, implying that their density must be 
< 7 x 10~ u pc~ 3 at the 95% confidence level. Comparing this limit to the density 
of local spheroid giants ~ 2 x 10~ 7 pc~ 3 (Morrison, 1993; Flynn & Fuchs 1994) 
we constrain the ratio of the densities of Local Group to galactic spheroid stars 
to be < 1/3000, about an order of magnitude lower than the limits obtained by 
Richstone et al. (1992). 
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